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ABSTRACT 


The  classical  equations  of  detonation  in  gaseous 
mixtures  are  presented  In  generalized  form  involving  only 
dimeni ionless  quantities.  This  report  sets  forth  a 
comprehensive  procedure  for  the  calculation  of  detonatlve 
Mach  numbers  and  associated  parameters  for  any  combustible 
gaseous  mixture. 

The  use  of  the  method  of  calculation  developed  is 
illustrated  by  means  of  a  sample  calculation.  For  this 
purpose,  calculations  are  presented  for  determination  of 
detonation  parameters  for  a  stoichiometric  mixture  of 
hydrogen  and  oxygen  at  an  initial  temperature  and  pressure 
of  18<>C  and  atmospheric  pressure,  respectively. 

The  calculation  method  developed  is  readily  adaptable 
to  computations  by  automatic  computing  machines  as  'well  as 
by  manually-operated  machines. 

Graphs  of  basic  data  are  presented  for  use  in  the 
sample  calculations,  and  it  is  believed  that  these  data 
will  find  practical  use  in  other  applications  as  well. 
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ments  involved  in  the  formation  of  one 
gram  mole  of  compound  "i",  relative 
to  the  standard  states  of  the  elements 
at  tL 

Jh^uj  Sum  of  the  enthalpies  (at  t^)  of  the  ele¬ 

ments  involved  in  the  formation  of  one 
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2 


M 


1 


Hv(AV  ^w(B)’  ^(CV 
^(D)-  ^w(E) 

Mw 

1 

Mw2 

Mw(i) 

Mw(j) 


Mw(U)  *  Mw(V)  *  Mw(W)’ 
Mw(X)*  Mw(Y)>  Mw(Z) 
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Mach  number;  Me  — 
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Detonative  Mach  number,  i.  e.  ,  the 
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with  respect  to  the  initial  mixture 


Mach  number  of  the  detonation  wave 
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Pseudo- Mach  number,  i.  e. ,  a  very 
close  approximation  to  the  detonative 
Mach  number,  M^;  used  for  conven¬ 


ience  only  and  defined  by 


2 

Mx  + 


(W* 


M, 


Molecular  weights  of  components  A, 

B,  C,D,E,  respectively,  in  the  ini¬ 
tial  mixture 

Average  molecular  weight  of  the  initial 
mixture 

Average  molecular  weight  of  the  final 
mixture 
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Molecular  weight  of  any  component 
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Number  of  moles 
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dimensionless 
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n. 


_•  _i  _i  _i  _i 

nA’VnC’VnE 
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Total  number  of  mole*  in  initial 
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Total  number  of  mole*  in  final 
mixture 

Number  of  mole*  of  any  component 
"i"  in  the  initial  mixture 

Number  of  moles  of  any  component 
"j"  in  the  final  mixture 

Number  of  mole*  of  components  A, 
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tial  mixture 

Number  of  mole*  of  components  U, 
V,W  X,  Y,  Z,  respectively,  in  the 
final  mixture 

Number  of  atoms  in  one  molecule  of 
a  given  component  of  a  mixture  (al¬ 
ways  an  integral  value  or  zero) 

Number  of  atoms  of  element  F,  G.  . . , 
respectively,  in  one  molecule  of  any 
component  Mi"  in  the  initial  mixture 
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Number  of  atoms  of  element  F,  G, . . . 
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Pressure 

Total  pressure  for  the  initial  mixture 
at  the  initial  conditions  (State  1) 

Total  pressure  for  the  final  mixture 
at  the  final  conditions  (State  2) 
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atm  os 
atmos 

atmos 
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SYMBOL 

Pi 


Pj 


DEFINITION 

Partial  prenure  of  any  component 
"i"  in  the  initial  mixture  at  the  ini¬ 
tial  conditions  (State  1) 

Partial  pressure  of  any  component 
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Pj  Total  pressure  of  the  initial  mixture 
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Total  pressure  of  the  final  mixture 
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ture  at  the  initial  conditions  (State  1) 


(Rj» 
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A  reference  temperature  pertaining 
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L  INTRODUCTION 


The  literature  on  the  subiect  of  detonation  in  gaseous  mixtures  dates  back  more 
than  seventy  years  11 1 -2-°a 2L HjAL.24 /  The  majority  o£  thia  iiterature  hag 

emphasised  the  prevention  of  detonation  rather  than  its  utilization.  In  connection  with 
combustion  in  propulsive  devices,  the  study  of  detonation  prevention  is  important  be¬ 
cause  of  the  destructive  effects  which  may  be  associated  with  uncontrolled  detonative 
combustion.  On  the  other  hand,  useful  application  of  detonative  type  combustion  in 
propulsive  devices  appears  to  warrant  consideration. 

In  the  normal  combustion  processes  occuring  in  the  combustion  chambers  of  pres¬ 
ent-day  combustion  devices,  the  flame  speeds  are  on  the  order  of  1  to  2  feet  per  second, 
whereas,  rates  of  propagation  for  stable  detonative  type  combustion  are  as  high  as 
3000  to  12,000  feet  per  second  for  various  gaseous  mixtures. 

In  order  to  be  able  to  reasonably  evaluate  the  possible  advantages  of  detonative 
type  combustion  in  propulsion  applications,  it  is  necessary  to  be  able  to  determine  the¬ 
oretically  the  values  of  all  detonation  parameters  for  various  gaseous  mixtures  con¬ 
sidered,  under  specified  initial  conditions.  Determination  of  these  detonation  para¬ 
meters  can  be  accomplished  experimentally;  however,  the  necessary  equipment  is 
highly  specialized  and  the  techniques  involved  are  not  sufficiently  developed  as  yet  to 
permit  the  obtaining  of  values  for  these  parameters  readily  under  all  desired  initial 
conditions. 

Methods  of  calculating  detonation  parameters  utilizing  directly  the  classical  deton¬ 
ation  equations  have  been  reported  by  Zeldovich  and  Ratnerii/  and  Edse  15/. 

In  performing  calculations  of  detonation  parameters  for  gaseous  mixtures,  it  was 
found  to  be  highly  desirable  to  rearrange  the  classical  (Chapman-Jouguet)  detonation 
equations  into  generalized  dimensionless  form.  In  this  report,  the  velocity  of  deton¬ 
ation  is  expressed  in  the  form  of  the  Mach  number  of  the  detonation  wave  relative  to  the 
unburned  gases;  the  heat  of  reaction  is  expressed  in  the  form  of  a  dimensionless  energy 
release  function.  The  other  detonation  parameters  are  all  expressed  as  dimensionless 
ratios  of  the  quantities  before  and  after  the  detonation  wave,  i.  e.  ,  pressure  ratio,  tem¬ 
perature  ratio,  density  ratio  and  molecular  weight  ratio.  Several  other  dimensionless 
quantities  are  introduced  to  simplify  further  the  graphical  presentation  of  the  generalized 
data  and  to  facilitate  calculations. 

A  method  is  presented  in  this  report  for  the  calculation  of  detonation  velocities 
and  other  detonation  parameters  for  gaseous  mixtures.  The  use  of  the  method  of  cal¬ 
culation  developed  is  illustrated  by  means  of  a  sample  calculation.  For  this  purpose, 
calculations  are  presented  for  determination  of  the  detonation  parameters  for  a  stoich¬ 
iometric  mixture  of  hydrogen  and  oxygen  at  an  initial  temperature  and  pressure  of 
18°C  and  atmospheric  pressure,  respectively. 

ROTS  -  This  Technical  Note  was  released  by  the  author  for  publication  in  Dec  saber  1954. 
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I.  INTRODUCTION 


The  literature  on  the  subject  of  detonation  in  gaseous  mixtures  dates  back  more 
17,18,19, . 


than  seventy  years  16. 


i  — — • - /  The  majority  of  this  literature  has 

emphasized  the  prevention  of  detonation  rather  than  its  utilization.  In  connection  with 
combustion  in  propulsive  devices,  the  study  of  detonation  prevention  is  important  be¬ 
cause  of  the  destructive  effects  which  may  be  associated  with  uncontrolled  detonative 
combustion.  On  the  other  hand,  useful  application  of  detonative  type  combustion  in 
propulsive  devices  appears  to  warrant  consideration. 


In  the  normal  combustion  processes  occuring  in  the  combustion  chambers  of  pres¬ 
ent-day  combustion  devices,  the  flame  speeds  are  on  the  order  of  1  to  2  feet  per  second, 
whereas,  rates  of  propagation  for  stable  detonative  type  combustion  are  as  high  as 
3000  to  12,000  feet  per  second  for  various  gaseous  mixtures. 

In  order  to  be  able  to  reasonably  evaluate  the  possible  advantages  of  detonative 
type  combustion  in  propulsion  applications,  it  is  necessary  to  be  able  to  determine  the¬ 
oretically  the  values  of  all  detonation  parameters  for  various  gaseous  mixtures  con¬ 
sidered,  under  specified  initial  conditions.  Determination  of  these  detonation  para¬ 
meters  can  be  accomplished  experimentally;  however,  the  necessary  equipment  is 
highly  specialized  and  the  techniques  involved  are  not  sufficiently  developed  as  yet  to 
permit  the  obtaining  of  values  for  these  parameters  readily  under  all  desired  initial 
conditions. 


Methods  of  calculating  detonation  parameters  utilizing  directly  the  classical  deton¬ 
ation  equations  have  been  reported  by  Zeldovich  and  Ratner— ^  and  Edse  15/. 

In  performing  calculations  of  detonation  parameters  for  gaseous  mixtures,  it  was 
found  to  be  highly  desirable  to  rearrange  the  classical  (Chapman- Jouguet)  detonation 
equations  into  generalized  dimensionless  form.  In  this  report,  the  velocity  of  deton¬ 
ation  is  expressed  in  the  form  of  the  Mach  number  of  the  detonation  wave  relative  to  the 
unburned  gases;  the  heat  of  reaction  is  expressed  in  the  form  of  a  dimensionless  energy 
release  function.  The  other  detonation  parameters  are  all  expressed  as  dimensionless 
ratios  of  the  quantities  before  and  after  the  detonation  wave,  i.  e.  ,  pressure  ratio,  tem¬ 
perature  ratio,  density  ratio  and  molecular  weight  ratio.  Several  other  dimensionless 
quantities  are  introduced  to  simplify  further  the  graphical  presentation  of  the  generalized 
data  and  to  facilitate  calculations. 

A  method  is  presented  in  this  report  for  the  calculation  of  detonation  velocities 
and  other  detonation  parameters  for  gaseous  mixtures.  The  use  of  the  method  of  cal¬ 
culation  developed  is  illustrated  by  means  of  a  sample  calculation.  For  this  purpose, 
calculations  are  presented  for  determination  of  the  detonation  parameters  for  a  stoich¬ 
iometric  mixture  of  hydrogen  and  oxygen  at  an  initial  temperature  and  pressure  of 
18°C  and  atmospheric  pressure,  respectively. 

NOTE  -  Ibis  Technical  Note  was  released  by  tbs  author  for  publication  In  December  1954* 
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IL  DERIVATION  OF  THE  GENERALIZED  EQUATIONS  IN  NON-DIMENSIONAL  FORM 


Some  of  the  assumptions  utilized  in  this  section  are  introduced  along  with  the  pre¬ 
sentation  of  particular  equations.  Other  assumptions,  more  general  in  nature,  will  be 
noted  before  proceeding  with  the  actual  derivation  of  the  equations. 

It  is  assumed  that  the  ideal  gas  law  holds  for  the  reactants  and  products  at  the  ini¬ 
tial  and  final  conditions,  respectively.  This  is  equivalent  to  the  assumption  that  the 
actual  temperatures  and  pressures  of  the  initial  and  final  mixtures  in  all  cases  are  in 
a  region  where  the  use  of  the  ideal  gas  law  results  in  negligible  error.  For  the  final 
mixture,  this  is  a  reasonable  assumption  in  all  cases  in  view  of  the  extremely  high 
temperatures  occuring  in  detonation.  For  the  initial  mixture,  this  is  a  reasonable  as¬ 
sumption  if  the  initial  conditions  are  sufficiently  far  from  the  critical  conditions  of 
pressure  and  temperature  for  each  component  so  that  compressibility  is  a  negligible 
factor.  Thus,  in  all  cases  where  an  equation  of  state  has  been  required  for  describing 
the  states  of  the  initial  and  final  mixtures  at  the  initial  and  final  conditions,  respec¬ 
tively,  the  ideal  gas  law  has  been  used. 

The  law  of  additives,  which  is  valid  for  mixtures  of  ideal  gases,  is  used  to  calcu¬ 
late  the  values  of  molecular  weight  (Mw)  and  molal  heat  capacity  (Cp)  for  the  initial 
and  final  mixtures. 

It  is  assumed  that  chemical  equilibrium  is  attained  in  the  detonation  wave. 

The  following  equations  of  conservation  across  the  detonation  front  are  applicable: 

Conservation  of  mass: 


p,  yi  *  P,yz 


(1) 


Conservation  of  momentum  (considering  only  normal  entry  into  the  front): 

(2) 


P,  +  pV,2*  P,  +  p  v,2 
1  r,  1  2  2 


Conservation  of  energy: 

V 2  p 

1  P1  T2 
Jel  +  T"  ■+  X  +  Jh  *  Je2  +  2  +  pt 


V,  p2 


(3) 


The  velocity  of  sound  in  the  initial  and  final  mixtures,  respectively,  may  be  ex¬ 
pressed  in  the  form 


■»*  •  U 


(4) 

(5) 
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The  assumption  is  made  here  that  the  adiabatic  law  holds  for  infinitesimal  changes  in 
the  initial  mixture  at  the  initial  conditions  and  in  the  final  mixture  at  the  final  conditions. 
This  is  merely  the  usual  assumption  made  in  deriving  the  expression  for  the  speed  of 
sound. 

Since  M  ■  V.  / a  and  M  *  V,/a  ,  then  with  equations  (4)  and  (5), 

1  A  1  2  *  C 

,  >  ->  X  M,2P, 

vi  *  MiS2  s  ~p  1  •  (6) 

*1 

v  2  w  2  2  T»  M22p2 

V  2.  *  ^2  *2  *  P  (^) 

rt 

By  definition,  the  energy  release  function  is 


From  equations  (8)  and  (4), 

Jh  «  Ba  2  .  B  •  (9) 

1  P 1 

From  the  conservation  of  momentum,  equation  (2)  and  equations  (6)  and  (7),  the 
pressure  ratio  across  the  detonation  front  may  be  expressed  as 


P2  „  ^1  ^1*  +  | 

7,  «.*  +  ' 


From  the  conservation  of  mass  equation  (1)  and  equations  (6)  and  (7), 


*  *  @  (t)  (i) 


and  from  equations  (10)  and  (11),  the  density  ratio  may  be  expressed  in  the  form 


A  t  y  M,*  \  /Y,  M,«  +  A 

^  *  \  T.M,**!/  I  T,  M,*  / 
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From  equations  (12)  and  (10), 


(Yi  M*  +l)* 


P»  P\  _  r«  Mi* 


p.  Pt  or,  Mf*  +o#  r.M,* 

The  equation  of  state  for  ideal  gases  may  be  expressed  in  the  following  form: 


(13) 


P  JRT 

~F*  ’ 


(14) 


Applying  equation  (14)  to  the  state  conditions  before  and  after  the  detonation  wave,  des¬ 
ignated  by  subscripts  1  and  2,  respectively,  the  following  relationship  is  obtained: 


P2  P\  T2Mw^ 
P1  Pt  *  Tl^w, 


(15) 


The  dimensionless  quantity — — - is  referred  to  hereafter  as  "the  temperature-mole* 

t1mw2 

cular  weight  ratio."  From  equations  (13)  and  (15),  the  following  expression  for  the 
temperature-molecular  weight  ratio  is  obtained: 


T,  M««  .  (Y,  M,f-f  l)* 

Ti  Mw,  (If,  Y,  M,» 


T,  Mw,  . 


(16) 


The  mass  and  momentum  conservation  equations  ((1)  and  (2))  have  been  utilized  to 


obtain  non-dimensional  expressions  for  -, 


P’  Pt 


1  Pt 


and 


T2Mw 


TiMw 


L  (equations  (10),  (12) 


and  (16)). 


In  order  to  similarly  transform  the  conservation  of  energy  equation  (3),  consideration 
of  the  internal  energy  term  is  required. 


The  first  law  of  thermodynamics  can  be  expressed  as 

dq  •  de  +  — dw  . 

el 


(17) 


For  the  case  where  the  only  force  acting  on  the  system  is  a  fluid  pressure,  dw  *  Pdv, 
and  equation  (17)  becomes 

dq  ■  de  +  JLPdv  .  (18) 

J 
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Differentiating  equation  (18)  with  respect  to  T,  holding  v  constant,  gives  the  following 
relationship: 


By  definition,  the  specific  heat  at  constant  volume  is 


and  it  follows  from  equation  (19)  that 


For  ideal  gases,  cv  is  independent  of  pressure  (or  volume)  and  therefore  is  a  function 
of  temperature  only.  Hence,  the  following  expression  may  be  obtained  from  equation 
(21): 


de 

(22) 

c  a 

V 

"dT"’ 

or 

de  a 

cvdT , 

(23) 

or 

Ae  a 

f  CydT  • 

(24) 

For  the  internal  energy  per  unit  mass  before  and  after  passage  through  the  detonation 
wave,  respectively,  equation  (24)  becomes 


el  ■  (^initial  mixture 


£ 


(cv).  ...  .  .  .  dT,  (25) 

*  vTlmtxal  mixture  '  ' 


*2  ~  ^e°Vinal  mixture 


£■ 


vVinal  mixture  dT> 


where  subscripts  1  and  2  refer  to  the  conditions  before  and  after  the  detonation  wave , 
respectively,  and  subscript  "o"  refers  to  the  reference  temperature  TQ.  Now,  intro¬ 
ducing  ^  as  a  mathematical  convenience,  equations  (25)  and  (26)  become 
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*1  "  *®o>initial  mixture  *  cv1#  0*T1  “  To*'  ^ 

*2  '  *e°Vinal  mixture  *  JTZ  '  To^  t28* 


where 


rl,o 


the  mean  specific  heat  at  constant  volume  for  the  initial  mixture 
between  state  1  and  the  reference  state  (T0), 


rZ,  o 


the  mean  specific  heat  at  constant  volume  for  the  final  mixture 
between  state  2  and  the  reference  state  (T0). 


The  reference  temperature  T0  will  be  taken  to  be  absolute  zero.  In  order  to  facilitate 
simplification  of  the  generalized  detonation  equations,  two  assumptions  will  be  made 
which  the  authors  do  not  feel  will  appreciably  affect  the  accuracy  of  the  final  general¬ 
ized  equations.  It  is  assumed  that  the  internal  energy  of  the  initial  mixture  at  the  ref¬ 
erence  temperature  (TQ  *  0®K)  is  equal  to  the  internal  energy  of  the  final  mixture  at 
the  reference  temperature  (TQ  *  O'K),  i.  e.  , 


^e°^initial  mixture  *  ^6°^final  mixture. 


The  mean  specific  heats,  7-  and  cv  ,  will  be  replaced  by  the  point  properties 

l.o  2,  o 

cv  and  cy  ,  respectively.  This  latter  substitution  or  assumption  will  permit  expres¬ 
sion  of  the^eneralized  detonation  equations  in  terms  only  of  the  point  properties  of  the 
initial  and  final  mixtures,  respectively. 

Based  on  the  foregoing  statements,  equations  (27)  and  (28)  may  be  combined  to 
yield  the  following  expression  for  the  difference  in  internal  energies  across  the  deton¬ 
ation  wave: 


'2  *  'l  *  Cv2T2  '  Vj  ’ 


(29) 


Utilizing  equation  (29),  the  conservation  of  energy  equation  (3)  becomes 


pi  vr 


P2  V2 

J°v2T2 

rz 


For  the  case  of  ideal  gases,  from  equation  (14), 


Mwi  P1  Mw2 

T1  *  11“  ’  ~p  and  T2  =  nr  ’ 


1 


(30) 


(31) 
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Using  aquations  (31),  the  energy  equation  (30)  may  be  written 


Since  cv  =  Cv^,  and  Cy^M^  *  Cy^,  and  -*nce»  for  ideal  gases, 
R  =  Cpi  -  CVi  =  Cp^  -  Cv^,  e<luat^on  (32)  becomes 

( JiL-V-i  ♦  !l*  ♦  \h  ♦  !l  . 

\cp,-SM  2  VVcvJ/>,  2 


(32) 


(33) 


and,  since  ^«Cpj/CVj,  and  ^*Cp^/Cv^,  the  following  equation  is  obtained  from 
eqviation  (33) : 


II. 

T,h 


Jh 


£  i 
T,-l  4  2 


(34) 


Equation  (34)  is  the  approximate  energy  equation  utilized  by  D.  G.  Samaras  in  his  paper 
on  endothermic  and  exothermic  discontinuities  (1 948)  LL/ and  involves  approximations  be- 
lieved  to  be  negligible,  as  discussed  above. 


From  equations  (6),  (7),  (9)  and  (34),  the  following  equation  is  obtained: 


JQ. 

Tr« 


/  Tfc  .  Yt  p»  Pi  , 

\Tg-i  2  /  p,  Pt 


and,  from  equations  (35)  and  (13), 


d  .  (  +  T&  .  IS  Mt* 

VT.-I  2  1 


(35) 


(36) 


This  is  a  form  of  the  Hugoniot  equation,  in  non-dimensional  form,  based  on  the  energy 
equation  (34). 


The  Chapman- Jouguet  hypothesis  that  the  product  gases  (final  mixture)  leave  the 
detonation  front  at  a  velocity  equal  to  the  local  speed  of  sound  in  those  gases  may  be 
expressed  in  the  form 

M2  =  1.0  .  (37) 
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With  this  substitution,  equation  (36)  reduces  to 


B  * 


/  tt«  .  ii*\  y« 

\r,-i  *  /  (w 


which,  upon  rearrangement,  gives  the  following  simplified  expression: 


(38) 


This  is  a  non-dimensional  form  of  the  equation  commonly  referred  to  as  the  Chapman- 
Jouguet  equation,  based  on  the  energy  equation  (34). 

Applying  the  Chapman- Jouguet  hypothesis  (equation  (37))  to  equations  $0),  (12)  and 
(16)),  expressions  for  the  pressure  ratio,  density  ratio  and  temperature-molecular 
weight  ratio  across  the  detonation  wave,  in  non-dimensional  form,  may  be  obtained 
as  follows: 


<*.”•♦'>  ■ 


(40) 


fi, 


* 


Ta  + 1  /  Yi  \ 

Y*  vt.  m,8+i/ 


(41) 


t,  mw,  _  Jz _  nr,  M(g  +o* 

T,  MW4  '  (T*+l)2  '  Yt  M,2 

The  above  three  equations,  together  with  equation  (39),  are  the  non-dimensional  gener¬ 
alized  detonation  equations  which  form  the  basis  for  calculations  and  graphical  presen¬ 
tations  in  this  report. 


Equations  (40),  (41)  and  (42)  are  presented  graphically  in  Figures  6,  7  and  8,  re¬ 
spectively.  In  each  of  these  three  graphs,  Mj2  has  been  made  the  abscissa.  Figure 
9  is  a  plot  of  Mj2  vs  Y,  Mj2  with  lines  of  constant  V|  and  is  presented  merely  to  show 
the  overall  range  of  the  term  X  M,2  encountered  in  the  detonation  of  gaseous  mixtures. 


T2Mw, 


j  U  W 1 

The  interrelationship  between  the  variables  P^/Pj,  P^/PK  an<*  " (t,ie  ^e_ 
pendent  variables  shown  in  equations  (40),  (41)  and  (42),  respectively)  is  presented  in 
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graphical  form  in  Figure  10.  It  will  be  noted  that  this  graph  ie  limply  a  plot  of  equa¬ 
tion  (15),  which  ii  the  equation  of  state  for  ideal  gases. 


IIL  REARRANGEMENT  OF  THE  CHAPMAN- JOUGUET  EQUATION 

The  Chapman- Jouguet  equation,  (39),  may  be  rearranged  and  expressed  in  the  fol¬ 
lowing  manner  in  order  to  facilitate  the  graphical  presentation  and  to  simplify  the  com¬ 
putational  procedure: 


M*  + 


(■*/  »*■ 2  i 


V-Y. 


nc-0 


(43) 


Equation  (43)  may  be  simplified  by  replacing  the  terms  on  the  left  hand  side  of  the  equa¬ 
tion  by  a  single  symbol,  M^,  and  replacing  the  two  terms  on  the  right  hand  side  of 

the  equation  by  simple  functions  of  the  dimensionless  parameters  £  and  P  ,  respec¬ 
tively.  The  term  M,  is  introduced  at  this  point  and  is  given  the  name  "pseudo-Mach 
number."  The  parameter  E  is  hereafter  referred  to  as  "the  dimensionless  parameter 
epsilon,  "  and  the  dimensionless  parameter  F*  is  hereafter  referred  to  as  "the  gamma 
function."  Thus,  equation  (43)  may  be  expressed  in  either  of  the  following  two  forms: 


M|*  -  2  (  r.'-O  B  +  2 


Y,  (V) 


(44) 


or 


M,2 


••  2  E  +  2f 


(45) 


where 


M,*  s 


M  2  + 


( _L 
'V  M,2 


(46) 


E  -  ( Ye*-l)  B 


(47) 
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(48) 


V-y 

Y,nr,-i) 

It  ihould  be  noted  that  equations  (43),  (44)  and  (45)  are  merely  rearrangements  of  the 
non-dimensional  form  of  the  Chapman-Jouguet  equation,  based  on  the  energy  equation 
(34). 

Equations  (47),  (48),  (45)  and  (46)  are  presented  graphically  in  Figures  1, 

2,3  and  4  respectively.  These  figures  present  the  non-dimensional  Chapman-  Jouget 
equation  in  graphical  form.  Figure  5  is  a  plot  of  V*.  vs  Yg  with  lines  of  constant 
and  is  presented  merely  to  show  the  overall  range  of  encountered  in  the  deton¬ 

ation  of  gaseous  mixtures. 


IV.  BASIC  DATA  REQUIRED 

A.  Data  on  heats  of  formation  for  the  reactants  at  the  initial  temperature  and  pres¬ 
sure,  or  these  data  at  a  standard  temperature  and  pressure  in  conjunction  with 
data  on  the  change  Cp  with  temperature  between  the  standard  conditions  and  the 
initial  conditions. 

B.  Data  on  molal  heat  capacity,  Cp^  ,  for  each  component  of  the  initial  mixture  at 
the  initial  conditions- 

C.  Data  on  molal  heat  capacity,  Cp(j),  for  each  component  of  the  final  mixture  at  the 
final  conditions  in  the  detonation  wave. 

D.  Data  on  equilibrium  constants  for  dissociation  reactions  involving  the  components 
of  the  final  mixture  at  the  final  conditions  in  the  detonation  wave. 


V.  DISCUSSION  OF  THE  CALCULATION  PROCEDURE 

The  calculation  procedure  based  on  the  generalized  detonation  equations  derived 
in  this  report  is  presented  in  detail  in  Appendix  1.  Immediately  preceding  this  pro¬ 
cedure  a  complete  outline  of  the  steps  involved  is  given  in  order  to  facilitate  its  use. 

Indications  are  that  previous  methods  for  calculating  detonation  parameters  have 
utilized  directly  the  classical  Chapman-Jouguet  equations  involving  dimensional  quan¬ 
tities,  whereas  in  this  report  the  calculation  procedure  utilizes  gener¬ 

alized  equations  and  non-dimensional  parameters  throughout. 
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This  method  of  calculation  involves  a  process  of  successive  approximations  or  a 
trial- and- error  procedure.  This  is  necessitated  by  the  fact  that  a  direct  solution  of 
the  generalized  equations  is  still  impractical,  if  not  impossible.  The  fact  that  disso¬ 
ciation  of  the  product  gases  is  taken  into  account  in  the  calculations  results  in  a  triple 
trial -and-err or  procedure.  This  procedure  involves  the  assignment  of  arbitrary  values 
to  final  pressure,  temperature  and  composition  and  then,  by  a  method  of  successive 
approximations,  the  values  of  these  quantities  are  adjusted  so  that  the  calculated  values 
approach  the  assumed  values.  Dimensionless  quantities  and  graphical  presentation  of 
the  generalized  data,  as  described  above,  are  used  to  simplify  and  streamline  consid¬ 
erably  the  calculation  procedure.  The  method  presented  provides  for  a  rapid  conver¬ 
gence  of  the  trial-and- error  calculations  by  simple  graphical  techniques.  The  calcu¬ 
lation  procedure  also  provides  for  continual  self-checking  of  results  at  intermediate 
points  throughout  the  entire  calculation,  i.  e.  ,  errors  can  be  detected  and  discarded 
with  little  loss  of  time. 

Considerable  effort  has  been  expended  in  the  preparation  of  this  method  to  insure 
that  the  information  contained  therein  is  clearly  understandable  to  any  research  scien¬ 
tist  who  might  have  occasion  to  be  working  in  this  field  (aeronautical,  mechanical, 
chemical,  electrical  or  propulsion  engineer,  or  chemist).  The  calculation  procedure, 
charts  and  graphs  are  presented  so  that  even  those  not  experienced  in  the  detonation 
field  can  readily  utilize  the  information  presented  to  obtain  accurate  values  of  the  re¬ 
quired  detonation  parameters  by  routine  manual  calculation  with  a  minimum  of  time, 
ingenuity,  instruction  and/or  supervision.  Where  extensive  data  at  different  initial 
conditions  and  for  various  gaseous  mixtures  are  desired,  this  method  is  also  very 
readily  adaptable  to  calculations  by  automatic  computing  machines. 


VI.  DISCUSSION  OF  THE  SAMPLE  CALCULATION 

A  sample  calculation  utilizing  the  calculation  procedure  presented  in  this  report 
is  given  in  Appendix  III.  Calculations  are  presented  for  the  determination  of  the  de¬ 
tonation  parameters  for  a  stoichiometric  mixture  of  hydrogen  and  oxygen  at  an  initial 
temperature  and  pressure  of  18* C  and  atmospheric  pressure,  respectively. 

Results  obtained  in  these  calculations  (taking  into  consideration  atomic  oxygen  in 
the  reaction  products  and  utilizing  consistent  data  throughout  for  specific  heats  and 
equilibrium  constants  in  the  range  of  the  final  temperature)  indicated  a  detonation  temp¬ 
erature  and  pressure  of  3536°K  and  17.  2  atmospheres,  respectively,  and  a  detonation 
velocity  of  8962  ft/ sec.  Results  obtained  by  this  calculation  procedure  neglecting  the 
dissociation  of  molecular  oxygen  but  still  utilizing  the  same  consistent  data  for  specific 
heat  and  equilibrium  constants  gave  values  of  the  detonation  temperature  and  pressure 
of  3593°K  and  17.  3  atmospheres,  respectively,  and  a  detonation  velocity  of  8991  ft/sec.* 

♦This  calculation  was  performed  by  the  authors,  but  detailed  calculations  are  not  pre¬ 
sented  in  this  report. 


WADC  TN  54-13 


11 


Calculated  data  obtained  from  literature  sources  indicate  a  detonation  temperature 
and  pressure  of  3583*K  and  18.05  atmospheres,  respectively,  and  a  detonation  velo¬ 
city  of  9203  ft/sec,  these  values  having  been  calculated  neglecting  the  dissociation 
of  molecular  oxygen^,  13/.  The  difference  between  the  calculated  data  presented  in 
this  report  and  that  reported  in  the  literature  is  probably  due  to  (a)  the  use  of  slightly 
different  specific  heat  and  equilibrium  constant  data  in  calculating  the  detonation  par¬ 
ameters,  and  (b)  the  simplifying  assumptions  utilized  in  deriving  the  generalized 
equations,  which  apparently  did  not  affect  greatly  the  accuracy  of  the  final  results. 

Preliminary  experimental  data  obtained  by  the  authors  in  a  one-inch  diameter 
detonation  tube  showed  very  close  agreement  (approximately  1-1/2%)  with  the  data 
calculated  by  the  procedure  reported.  The  experimental  value  obtained  for  the  deton¬ 
ation  velocity  was  9100  ft/sec. 
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APPENDIX  I 


CALCULATION  PROCEDURE 


Detailed  Outline  of  Calculation  Steps 
Detailed  Calculation  Procedure 
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For  convenience  in  using  this  calculation  procedure,  a  detailed  outline  of  the  steps 
involved  is  presented  below. 

Step  1.  Proposed  Scheme  of  Reaction 

Step  2.  Setting-Up  of  the  Partial  Pressure  Equations  for  Use  in  the  Calculation 
of  the  Composition  of  the  Final  Mixture 

Step  3.  Solution  of  the  Simultaneous  Partial  Pressure  (pj)  Equations  in  General 
Terms 

Step  4.  Fixing  of  the  Initial  Composition,  n^1 

Step  5.  Particular  Solution  of  the  Partial  Pressure  (pj)  Equations  for  a  Specific 
Initial  Mixture  Composition  (n£) 

Step  6.  Fixing  of  the  Initial  Pressure,  p^ 

Step  7.  Fixing  of  the  Initial  Temperature,  t^ 

Step  8.  Calculation  of  the  Properties  of  the  Initial  Mixture 

Step  9.  Assumption  of  the  Final  Pressure,  p^ 

Step  10.  Assumption  of  the  Final  Temperature, 

Step  11.  Determination  of  Numerical  Values  of  the  Equilibrium  Constants  for  the 
Dissociation  Reactions 


Step  12.  Calculation  of  the  Partial  Pressures  (p.'s)  in  the  Final  Mixture 

J 

Step  13.  Determination  of  the  Total  Number  of  Moles  (n^)  of  the  Final  Mixture^ 
Step  14.  Calculation  of  the  Composition  (nj's)  of  the  Final  Mixture 
Step  15.  Calculation  of  the  Properties  of  the  Final  Mixture 
Step  16.  Calculation  of  the  Heat  of  Reaction,  A 


Step  17.  Calculation  of  the  Energy  Release  Function,  B 
Step  18.  Determination  of  the  Value  of  the  Gamma  Ratio, 


Step  19.  Calculation  of  the  Molecular  Weight  Ratio,  Mw 

Step  20.  Determination  of  the  Value  of  the  Dimensionless  Parameter 

Step  21.  Determination  of  the  Value  of  the  Gamma  Function,  P 
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Step  22.  Determination  of  the  Value  of  the  Paeudo-Mach  Number  (Squared),  Mj2 

Step  23.  Determination  of  the  Value  of  the  Detonative  Mach  Number 
(Squared),  M, 


Step  24.  Determination  of  the  Value  of  Mj2 

Step  25.  Determination  of  the  Value  of  the  Temperature-Molecular  Weight  Ratio, 

T£Mw^  /TiMw^ 

Step  26.  Calculation  of  the  Final  Temperature,  T^ 

Step  27.  Determination  of  the  Value  of  the  Pressure  Ratio,  p  /p 

£  i 

Step  28.  Calculation  of  the  Final  Pressure,  p 

£ 

Step  29.  Calculation  of  the  Detonative  Mach  Number,  M^ 

Step  30.  Calculation  of  the  Detonation  Velocity,  D 


Step  31.  Calculation  of  Additional  Data 

a.  Density  ratio,  P^f  P\ 

(or  Volume  Ratio,  Vj/v ^) 


b.  Initial  Density,  p  (or  — — ) 

c.  Final  Density,  p  (or  ) 

t  2 

The  detailed  calculation  procedure  is  presented  below  and  follows  the  outline  given 
directly  above. 


Step  1.  Proposed  Scheme  of  Reaction 

Set  down  initially  a  proposed  scheme  of  reaction,  including  all  pos¬ 
sible  or  probable  products: 

■'aA  +  '<BB+ - ►  nuU  +  »vV+...,  (49) 

ri ^  and  rt  being  known  quantities,  n^  and  n  being  initially  unknown  quan¬ 
tities.  The  number  of  possible  products  included  in  this  reaction  scheme 
will  depend  upon  experience  and  information  available  on  the  particular 
reaction  being  investigated. 

The  calculation  procedure  should  be  started  on  the  basis  of  the  above 
proposed  reaction  scheme.  After  the  first  calculation  of  composition, 
based  on  the  first  assumed  temperature  and  pressure,  it  may  be  apparent 
that  one  or  more  of  the  products  in  the  proposed  reaction  scheme  is 
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present  only  in  negligible  amount.  If  this  is  the  case,  then  it  may  be 
advisable  to  modify  the  above  proposed  scheme  of  reaction  to  eliminate 
the  negligible  constituents  before  proceeding  with  the  second  calculation 
of  composition. 


Step  2.  Setting-Up  of  the  Partial  Pressure  (pj)  Equations  for  use  in  the  Calcula-. 
tion  of  the  Composition  of  the  Final  mixture 


a.  Dissociation  Reactions  and  Equilibrium  Equations 


Set  down  dissociation  reactions  involving  each  of  the  products 
included  in  the  above  reaction  scheme,  along  with  the  corresponding 
equilibrium  equations  expressed  in  terms  of  the  partial  pressures 
of  the  constituents.  For  example,  if  ^O,  H2  and  are  three  of  the 

components  of  the  final  mixture,  then 


H  O  H  +  1/2  O 

6  2  “ 


and 


Kp  <h2o) 


equilibrium  constant. 


(50) 


(51) 


In  selecting  the  dissociation  reactions  to  be  used  in  these  calcu¬ 
lations,  care  must  be  excercised  to  insure  that  the  equilibrium  equa¬ 
tions  obtained  are  independent  equations.  Each  successive  dissocia¬ 
tion  reaction  written  must  involve  a  new  component  of  the  final  mix¬ 
ture  if  the  corresponding  equilibrium  equation  is  to  be  independent  of 
the  other  equilibrium  equations.  Another  factor  to  be  considered  in 
the  selection  of  the  dissociation  reactions  to  be  used  is  the  availabil¬ 
ity  of  equilibrium  constant  data  for  particular  reactions. 

b.  Additional  Equations  Required 

The  dissociation  reactions  have  provided  a  number  of  independent 
equations  involving  the  partial  pressure  of  each  component  of  the  final 
mixture.  In  order  to  solve  for  each  of  these  partial  pressures,  the 
total  number  of  independent  equations  available  must  equal  the  number 
of  components  in  the  final  mixture.  Independent  equations  required 
in  addition  to  the  equilibrium  equations  are  obtained  from  (1)  applica¬ 
tion  of  Dalton's  Law,  i.  e.  ,  summation  of  partial  pressures,  and  (2) 
application  of  the  law  of  conservation  of  mass  for  each  element  in¬ 
volved  in  the  proposed  scheme  of  reaction. 
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c.  Summation  of  Partial  Pressure*  in  the  Final  Mixture 

Set  down  the  equation  for  the  summation  of  the  partial  pressures 
of  the  components  of  the  final  mixture: 


PU  +  PV  +  PW  +-“  ** 


d.  Conservation  of  Mass  for  Each  Element 


(52) 


Set  down  the  mass -balance  equation  stating  the  conservation  of 
atomic  type  for  each  chemical  element  involved  in  the  proposed  reac¬ 
tion  scheme.  The  number  of  these  mass-balance  equations  is  equal 
to  the  number  of  chemical  elements  involved.  The  law  of  conserva¬ 
tion  of  mass  yields,  for  elements  F  and  G,  the  following  two  equations: 

For  element  F: 

if  ^  t  nig  Ng  +  • . .  —  ny  +  n^.  Ny  + . . . ,  (53) 

or 

X  ri/Ni  =  £  nj'Nj  .  (54) 

For  element  G: 

riA*NA  +  rig9NB  +. .  .  a  nu*^U  +  nV#NV  +*  •  •  *  (55) 

or 

X  rii*Ni  *  2nj®Nj  •  <56> 

i  j 

Divide  each  mass  balance  equation  by  n^  (value  initially  unknown) 
as  shown  below.  For  element  F,  the  equation  obtained  is 


n. 


i  2 


(57) 


and,  similarly,  for  element  G, 


19 


(58) 


WADC  TN  54-13 


Substitute  into  the  mass  balance  equation  for  each  element  the 
following  equality,  obtained  from  Dalton's  Law  of  Additive  Pressures 
and  the  ideal- gas  law: 


nj  pj 

■  ■  »  *5  • 

n2 


(59) 


This  yields  a  mass  balance  equation,  for  each  element,  in  terms  of 
the  partial  pressures  (Pj's)  of  components  of  the  final  mixture.  For 

element  F,  this  yields  the  equation 


(60) 


and,  for  element  G,  the  equation 


or,  rearranging  slightly,  for  element  F, 


and,  for  element  G, 


(61) 


(62) 


rit 


(63) 


As  stated  above,  the  number  of  these  mass  balance  equations 
(equations  (62)  and  (63))  is  equal  to  the  number  of  chemical  elements 
involved  in  the  proposed  reaction  scheme.  Any  two  of  the  mass  bal¬ 
ance  equations  may  be  combined  to  yield  an  independent  equation  in 
terms  of  the  partial  pressures  (pj's),  by  elimination  of  n^  and  p^.  For 
example,  for  elements  F  and  G,  dividing  equation  (62)  by  (63)  yields 

2  n.'  'n,  2  p  V 

1  -1—1  =  i  A  .  (64) 

?  n;  *N,  2  p.  *Nj 
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If  there  are  three  chemical  elements  involved  in  the  reaction 
scheme,  an  additional  independent  equation  in  terms  of  the  partial 
pressures  (pj's)  may  be  obtained  by  using  a  new  combination  of  mass 
balance  equations.  If  more  than  three  elements  are  involved,  addi¬ 
tional  independent  equations  in  terms  of  the  partial  pressures  may  be 
obtained  in  the  same  manner. 

Step  3.  Solution  of  the  Simultaneous  Partial  Pressure  (pj)  Equations  in  General 
Terms 

Step  2  has  provided  a  number  of  independent  partial  pressure  equa¬ 
tions  equal  to  the  number  of  components  in  the  final  mixture.  Solve 
these  equations  simultaneously  to  obtain  a  general  solution.  When  the 
number  of  products  of  reaction  is  small,  this  solution  of  the  simultaneous 
equations  may  be  accomplished  by  direct  algebraic  methods.  When  the 
number  of  reaction  products  is  large,  it  may  be  desirable  to  utilize  meth¬ 
ods  of  solution  such  as  those  reported  by  Edse.-S^  Huff  and  Calvert,  ^ 
Huff  and  Morrell,-^  Hershfelder , etc. 

Particular  solutions  for  specific  initial  mixtures  (rfj)  are  obtained 
from  the  general  solution  by  fixing  the  initial  mixture  composition  (Steps 
4  and  5). 

Step  4.  Fixing  of  the  Initial  Mixture  Composition,  rij 

Fix  the  initial  mixture  composition,  n^,  at  the  value  to  be  used 
throughout  the  remainder  of  the  calculations. 

Step  5.  Particular  Solution  of  the  Partial  Pressure  (pj)  Equations  for  a  Specific 
Initial  Mixture  Composition  (nfj) 

Substitute  the  numerical  value  for  the  initial  mixture  composition 
(nfj),  assigned  inStep  4,  into  the  general  solution  obtained  in  Step  3,  to 
obtain  a  particular  solution  of  the  partial  pressure  (pj)  equations  for  a 
specific  initial  mixture  composition. 

Step  6.  Fixing  of  the  Initial  Pressure,  p^ 

Fix  the  initial  pressure,  p^,  at  the  value  to  be  used  throughout  the 
remainder  of  the  calculations. 

Step  7.  Fixing  of  the  Initial  Temperature,  t^ 

Fix  the  initial  temperature,  t^,  at  the  value  to  be  used  throughout 
the  remainder  of  the  calculations. 

Step  8.  Calculation  of  the  Properties  of  the  Initial  Mixture 
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a.  Mole  Fraction,  xfj 


(65) 


x'  • 


nA  +  nB  +  *‘* 


or 


5  n'i 


b.  Molecular  Weight,  Mw 


1 


MWj  =  x'aMw(A)  +  x'b  ^w(B)  +  • 
=  2  x'i^w(i) 


(66) 


(67) 


c.  Molal  Heat  Capacity  at  Constant  Pressure,  C_ 


d. 


e. 


xA^p(A)  +  x'b^p(B) 


+  . . . 


•  2*ic'P<i>  • 

1 

Specific  Heat  Ratio  (or  Molal  Heat  Capacity  Ratio), 


Heat  of  Formation*,  A  H'f^ 


T, 


(68) 


(69) 


Calculate  the  heat  of  formation  of  each  component  of  the  ini¬ 
tial  mixture  at  the  initial  conditions. 

AH'f(i)  =  AH,^  -H'i+  2  H„(i)  .  (70) 

If  the  component  "i"  is  an  element  in  its  standard  state  of  aggre¬ 
gation  at  tj ,  then,  by  definition,  it  is  apparent  that  its  heat  of 

formation  AH'f^  is  equal  to  zero,  e.g. ,  j  and  A  H'^q  j 

2  2 

are  equal  to  zero  when  tj  is  in  a  region  where  hydrogen  and 

*  In  this  report,  the  heat  of  formation  is  taken  to  be  positive  when  heat  is  evolved.  If 
enthalpy  data  for  this  calculation  cannot  be  found  in  the  literature,  it  can  be  calculated 
readily  from  available  Cp  data. 
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Step  9. 


Step  10. 


Step  11. 


Step  12. 


Step  13. 


oxygen  are  normally  diatomic  gases. 

Assumption  of  the  Final  Pressure,  p.;s 

Assume  a  value  for  the  final  pressure,  p~,  . 

Assumption  of  the  Final  Temperature, 

Assume  a  value  for  the  final  temperature,  • 

Determination  of  the  Numerical  Values  of  the  Equilibrium  Constants  for 
the  Dissociation  Reactions 

Obtain  the  numerical  values  of  the  equilibrium  constant  for  each  of 
the  dissociation  reactions  considered  in  Step  2,  at  the  assumed  temp¬ 
erature  (T^),  either  from  available  data  in  the  literature  or  by  calcu¬ 
lation. 


Because  of  the  trial- and- error  nature  of  the  overall  calculation  pro¬ 
cedure,  it  is  advisable  to  make  a  plot  of  equilibrium  constant  versus 
temperature  for  each  of  the  dissociation  reactions  in  the  temperature 
range  in  which  *8  expected  to  be.  It  is  necessary  that  this  be  done 

at  the  very  beginning  of  the  calculation  procedure  in  order  to  insure 
that  consistent  data  are  used  and  thereby  facilitate  rapid  convergence 
of  the  trial- and- error  calculations. 


Calculation  of  the  Partial  Pressures  (pj's)  in  the  Final  Mixture 

Using  tho  squat ion(s)  obtained  in  Step  5  in  conjunction  with  the  equil¬ 
ibrium  equations  set  forth  in  Step  2,  par.  a,  determine  the  numerical 
value  of  the  partial  pressure  (pj)  of  each  constituent  of  the  final  mixture. 

Determination  of  the  Total  Number  of  Moles  ( n ^  of  the  Final  Mixture 

Set  down  the  overall  mass  balance  equation  for  the  proposed  scheme 
of  reaction: 

“'A^MA)  +  n'BMw(B)  +  •••  =  nUMw(U)  +  “V^MV)  +  •••  »  (71> 


or 


2  n  i  ^w(i)  *  2  nj  * 

i  j 

Substitute  in  equation  (72)  the  relationship 


(72) 


(73) 
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obtained  from  equation  (59);  then 


2  n'i  Mw(i)  =  T”  2  Pj  ^w(j)  ' 

i  y2  j 


and 


n'iMw(i) 
Pj  Mw(j) 


(75) 


This  gives  the  total  number  of  moles  of  final  mixture  corresponding  to 
the  number  of  moles  of  initial  mixture  shown  in  the  proposed  scheme  of 
reaction  (i.  e.  ,  n'^  +  n'g  +  ...  ). 

Step  14.  Calculation  of  the  Composition  (nj's) 

Again  using  the  relationship 


n. 


(73) 


obtained  from  equation  (59),  calculate  the  number  of  moleB  (n.)  of  each 
component  in  the  final  mixture. 

Step  15.  Calculation  of  the  Properties  of  the  Final  Mixture 


a.  Mole  Fraction,  x- 

J 


-  JlH.  - 

XU  “  n2  '  p2 


b.  Molecular  Weight,  Mw 

Cm 

Mw  *  xUMw(U)  +  xVMw(V)  +•••  * 

Cm 


2xjMw(j)  • 
j 


(76) 


(77) 


(78) 
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c.  Molal  Heat  Capacity  at  Constant  Pressure,  Cp^ 

Determine  the  value  of  C„  at  the  above  assumed  value  of  T, 

P2  2 

(Step  10),  using  values  of  Cp^j  obtained  from  the  literature. 

Cp2  *  *uCp(U)  +  xVCp(V)  +•••  =  2xicp(j)*  <79> 

J 

d.  Specific  Heat  Ratio  (or  Molal  Heat  Capacity  Ratio) , 


P2 


P2 


Cp2-R 


(80) 


If  Cp  data  cannot  be  found  in  the  literature,  it  can  be  readily  calculated 
from  available  enthalpy  data. 


e.  Heat  of  Formation*,  AH' 


f(j) 


Calculate  the  heat  of  formation  of  each  component  of  the  final  mix¬ 
ture  at  the  initial  conditions. 


.  AH,  (j)  H'j  +  2H'e(j) 


(81) 


Step  16.  Calculation  of  the  Heat  of  Reaction,  Energy  Release  per  Mole  of  Initial 
Mixture  (AH^) 

nU  AH’f(u,  +  ny  AH'jjyj  +...  -  nA  ^H'f(A)  "  nB  AH'f(B) 

AHR  *  nA  +  nfi  +...  * 

(82) 

or 

£  nj  AH',(j)  -  Sn'iAH',,,, 

AH„  *  — - ^ -  *  (83) 


Step  17.  Calculation  of  the  Energy  Release  Function,  B 

By  equations  (8),  (4)  and  (31)  and  by  definition  of  "h"  and  AH^, 


*  See  footnote  in  Step  8,  par.  e. 
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A  Hr 


„  _  Jh 

^  *  ttrt,  *  Trt;  • 


(84) 


Step  18.  Determination  of  the  Value  of  the  Gamma  Ratio, 

T 

Calculate  this  ratio  from  values  of  Kj  and  Tfjor  obtain  from  Figure  5. 


Step  19.  Calculation  of  the  Molecular  Weight  Ratio, 


mW2 

Mw, 


Calculate  this  ratio  from  values  of  Mw  and  1L,  . 

W1  Z 

Step  20.  Determination  of  the  Value  of  the  Dimensionless  Parameter  £ 

Calculate  £  from  values  of  K^and  B  or  obtain  from  Figure  1. 

E»  b 

Step  21.  Determination  of  the  Value  of  the  Gamma  Function,  P 

Calculate  P  from  values  of  2T,  and  Tg  or  obtain  from  Figure  2. 


(47) 


r  ■ 


(48) 


TT  2 


Y,(V) 

Step  22.  Determination  of  the  Value  of  Pseudo-Mach  Number  (Squared),  Mj‘ 
Calculate  Mj2  from  values  of  E  and  T  or  obtain  from  Figure  3. 

Mj2  =  2E  +  ,  (45) 


or 


x  v-l 

M..2(Ti-0B+2__. 


(44) 


Step  23.  Determination  of  the  Value  of  the  Detonative  Mach  Number  (Squared),  Mj2 
Obtain  the  value  of  M^  from  Figure  4. 

n uxf 


M,  *  M, 


(46) 


M, 
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Step  24.  Determination  of  the  Value  of  K|  M(* 

Calculate  the  value  of  Kj  M*  from  values  of  T(  and  M,*  or  obtain 
from  Figure  9. 

Step  25.  Determination  of  the  value  of  the  Temperature-Molecular  Weight 
T2Mw  j 


Ratio, 


Tl“w2 


T2Mw  .  -y  2  Y* 

Calculate  the  value  of  - — — = —  from  values  of  0.  M,  and  0» 

HP  a  a  •  *  * 


t1Mw2 

or  obtain  from  Figure  8.  From  equation  (42), 


T2  Mwi 


r. 


(Ta 


+  ,)*  (*•  "'*  +  Y,  M,*+2) 


(85) 


Step  26.  Calculation  of  the  Final  Temperature,  T 2 

-t-wj  i»w> 

Calculate  the  value  of  T-,  from  values  of  — -  ,  - =  and  T,. 

2  TlMw,  Mwj  1 


t2mWi 


T 


2 


T2MWj  MW£ 
T1Mw2 


(86) 


If  this  calculated  value  of  T2  does  not  check  the  assumed  T2  (Step  10), 
re-assume  T2  and  repeat  Steps  10  through  26.  In  order  to  facilitate  rapid 
convergence  between  calculated  and  assumed  values,  the  following  pro¬ 
cedure  may  be  used.  After  two  calculations  of  T2,  plot  calculated  T2 
versus  assumed  T2  and  draw  a  line  through  these  two  points,  intersecting 
the  45#  line.  The  value  of  T2  at  the  point  of  intersection  usually  is  a  very 
good  indication  of  the  correct  value  of  T2  corresponding  to  the  assumed 
p2<  This  value  of  T2  should  be  used  as  the  assumed  T2  for  third  trial.  A 
typical  plot  is  shown  in  Figure  11. 

P2 

Step  27.  Determination  of  the  Value  of  the  Pressure  Ratio,  — 

Pi 


Calculate  the  value  of  ~  using  the  values  of  Kg  and  "K|  Mj* 

Pi 


corres¬ 


ponding  to  the  final  checked  value  of  T2  obtained  in  the  trial-and- error 
calculations  of  Steps  10  through  26,  or  obtain  from  Figure  6.  From 


equation  (40), 


since  is  numerically  equal  to 


II 

Pi  ’ 
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£2 

Pi 


Step  28.  Calculation  of  the  Final  Pressure, 


(87) 


(88) 


If  thia  calculated  value  of  p^  does  not  check  the  assumed  p^  (Step  9),  re¬ 
assume  p^  and  repeat  Steps  9  through  28  until  the  calculated  and  assumed 
values  agree. 


The  graphical  procedure  outlined  under  Step  26  may  be  used  to  facil¬ 
itate  rapid  convergence  between  calculated  and  assumed  values  of  p  . 

Step  29.  Calculation  of  the  Oetonative  Mach  Number, 

Calculate  Mj  from  the  value  of  (Step  23)  corresponding  to  the 
final  checked  values  of  p^  and  T ^  (Steps  26  through  28). 

Step  30.  Calculation  of  the  Detonation  Velocity,  D 

Since  D  *  *  Mja^  then  with  equations  (4)  and  (14), 


D  =  Mj 


(89) 


Step  31.  Calculation  of  Additional  Data 

/>, 


1 


a.  Density  Ratio,  JJL  (or  Volume  Ratio, - ) 

P\  V2 

D  v  p 

Calculate  (or  — — )  from  the  final  checked  values  of  — — 

Py  v2  Pi 

T2MW| 

and — —  — - obtained  in  the  trial-and-error  calculations  of  Steps  26 

through  28,  or  obtain  from  Figure  10. 


b.  Initial  Density, 


(15) 


W  ADC  TN  54-13 


28 


By  equation  (14),  for  the  initial  mixture, 


c. 


MwjPj 

JRT, 


Final  Density,  p  (or  -JL  ) 

^  2 

p 

P9  r  JL  .  p 

'2  ft  n 


(90) 


(91) 
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APPENDIX  II 


GRAPHS  OF  BASIC  DATA-  - 
EQUILIBRIUM  CONSTANTS  AND 
MOLAL  HEAT  CAPACITIES 
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As  pointed  out  in  the  Calculation  Procedure  in  this  report,  it  is  desirable  that 
consistent  data  on  equilibrium  constants  and  molal  heat  capacities  be  available.  Since 
the  sample  calculation  presented  in  this  report  (Appendix  III)  is  for  a  stoichiometric 
hydrogen- oxygen  system,  the  necessary  equilibrium  constant  data  and  molal  heat  cap¬ 
acity  data  for  this  system  are  presented  in  this  appendix. 

Values  of  the  equilibrium  constants  (K  's)  as  a  function  of  the  absolute  tempera¬ 
ture  for  H2O,  H2  and  O2  are  presented  graphically  in  Figures  12  through  15,  respec¬ 
tively. 

Values  of  the  molal  heat  capacities  (Cp's)  as  a  function  of  absolute  temperature 
for  H2O,  H2.  O2  and  OH  are  presented  graphically  in  Figures  16  through  19,  respec¬ 
tively. 

The  range  of  temperature  in  all  these  figures  extends  from  2000°  to  5000* K,  and 
it  is  expected  that  this  range  will  cover  all  cases  to  be  encountered  in  the  detonation 
of  gaseous  mixtures. 

In  each  of  the  figures  described  above,  a  smooth  curve  was  drawn  through  the 
data  plotted  in  order  to  insure  that  consistent  data  would  be  available  for  the  trial-and- 
error  calculations. 

It  is  believed  that  these  graphs  of  basic  data  described  above  will  find  practical 
use  in  other  applications,  particularly  in  other  combustion  calculations  where  trial - 
and-procedures  are  involved  and  where  consistent  thermodynamic  data  are  essential. 
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APPENDIX  III 


SAMPLE  CALCULATION 
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CALCULATION  OF  DETONATIVE  MACH  NUMBER  AND  ASSOCIATED  PARAMETERS 

FOR  H2  -  Oz  MIXTURES 

The  step  numbers  used  below  are  the  same  as  those  given  in  the  calculation  pro¬ 
cedure  outlined  in  Appendix  I. 

Step  1.  Proposed  Scheme  of  Reaction 

n*AA  +  n'BB“*  nyU  +  nvV  +  n^W  +  n^X  +  nyY  +  n2Z  (92) 


or,  specifically,  for  H2  -  02  mixtures, 

nHtH*+  n0,  OST*'\a0f,MO  +  +  "08°*  +  "OH011  +  “hH  +  ”6°  • 

Step  2.  Setting-Up  of  the  Partial  Pressure  (pj)  Equations 
a.  Dissociation  Reactions 

(i)  h2o  5=±  h2  +  io2 


kP(M#0) 


PHf  *  V 


’to 


(2)  HzO 


P(HtO) 


jHz  +  OH 

^  ■  P.H 

V 


<3>  i»z 


H 


1/  .  n 

Kp  “  i 

*"«'  V 

(4)  ^o2  3=ato 

„  .  fo 

Kp-v  ’  !$' 

b.  Summation  of  the  Partial  Pressures  in  the  Final  Mixture 


(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

(100) 

(101) 


Pu  +  PV  +  PW+  PX  +  PY  +  PZ  *  I  Pj  *  Pt 

j 


(102) 
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or,  specifically,  for  Hg  -  O %  mixtures, 

PH«0  *  PM,  +  P0*  +  P0H  +  PH  +  P0  *  P* 

c.  Conservation  of  Mass  for  Each  Element 

Using  equation  (64), 

2  n!  FNt  2  p.  rNj 

?  pi  Ni 


(103) 


S  n\  #N 


(64) 


i 


and  taking  elements  F  and  G  to  be  O ^  and  H^,  respectively,  then 


>Z  PM«0  +  2  P0.  +  P0H  +  P0 
—  S  - = - = - 


'Hf 


^PHtO  +  ^PH|  ■*"  P0H  +  PH 


(104) 


or 


(2  -•)  v+ 2 vrt  p-,+ ifct  p-+(^-')  po«“2po,_po=0- <io5> 


Step  3. 


Solution  of  the  Simultaneous  Partial  Pressure  (pj)  Equations  in  General 
Terms 


Step  2  has  provided  the  following  six  independent  partial  pressure 
equations,  which  are  sufficient  to  determine  the  partial  pressures  of  each 
of  the  six  components  in  the  final  mixture. 


K 


<HS0) 


s 


(95) 


(HgO) 


(97) 


Pu 

Kp  *  - V-  *  (99) 

pT 
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Kn  *  ^  . 

<o*>  p.I 


(101) 


PMf0  +  PHg  +  POf  +  P0M  +  PH  +  P0  s  Pt  ’ 
(  n'°t  \  n0*  n°*  fnOm  \ 

Rearranging  equations  (95),  (97),  (99)  and  (101), 


(103) 


(105) 


iHgO) 

P 


_VY 

Ht  / 


(106) 


POH  5 


^P<HgO)  Ph2° 


P  1 


Pn=  Kn  Prt 

0  P(0*>  °* 


Kn  P  *  , 

P(Hg>  M* 

i  Kp(Og)  Kp(HgO)  PHgO 

■  z  . - . —  t 


(107) 


(108) 


(109) 


and  substituting  these  relationships  into  equations  (103)  and  (105),  re¬ 
spectively,  the  following  two  independent  quadratic  equations  in  terms 

of  two  variables,  n  and  p  ,  are  obtained: 

TIjO  Mg 


KP(MgO>Y  Jt  .  /KP(HgO)  ,  KP<0*)  ^PlMgO) 


Ph-O  + 


■*)***  (WlPM*i"P‘)S0,  (110, 


{*-*&.)'  p*  +  [“ 8  V  J  pHJo 
\  pm  /  \_  pH,  '"«*  1  P„*  n"«  J 

-^(2V  K%>p*)s°- 


(111) 
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Step  4.  Initial  Composition 


For  the  purpose  of  this  sample  calculation,  a  stoichiometric  mix¬ 
ture  of  hydrogen  and  oxygen  will  be  considered,  the  initial  composition 
being  given  by 


where 


n«,H.  +  . 


and  n 


o*  -  2  » 


or 


^  +  ^02  . 


(112) 


Step  5.  Particular  Solution  of  the  Partial  Pressure  (pj)  Equations 


From  the  values  for  initial  composition  given  in  Step  4, 


n 


°*  _ 


Hg 


(113) 


Substitution  of  the  numerical  value  of  this  ratio  into  the  general  partial 
pressure  equations  (Equations  (110)  and  (111))  obtained  in  Step  3  yields 
the  following  particular  solution  for  a  stoichiometric  mixture  of  hydrogen 
and  oxygen: 


and 


(IS-ulY  P.  . 
\  %)  M*° 


KP(0«)  *P(HtO>  +  KP(HgO) 


H« 


"K+iKP, 


j 


HgO 


P  I  ): 
'«>  H*  / 


o. 


(114) 


(115) 
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Step  6.  Initial  Pressure 


For  the  purpose  of  this  sample  calculation,  the  initial  pressure  is 
fixed  at  1  atmos. ,  i.  e. , 


and 


Pj  =1  atmos. 


P  =  1.0133  x  10^  dynes/cm^. 
1 


Step  7.  Initial  temperature 


For  the  purpose  of  this  sample  calculation,  the  initial  temperature 
is  fixed  at  18*C,  i.  e.  , 


and 


tj  ■  18°C. 
T  ■  291°K. 


Step  8.  Properties  of  the  Initial  Mixture 

a.  Mole  Fractions 


1 

.  1  . 

Vn0t 

'  •+?  ' 

Z  5 

n'0l 

.  T 

1 

nH,+  nO, 

,  +  i 

*  3  » 

^  x, .  ><,,  +  K, *  '  • 

b.  Molecular  Weight 

MW,  S  *Ht  MW(h#v  +  *0*  MW<0f) 
*  (J)  (2.0160)  +  (3)  (32.000) 


•  12.0107. 

c.  Molal  Heat  Capacity 
At  tj  ■  18#C.  , 

c'n  ■  6.  88-^  cal/(gm  mole)(#K); 
K(Hf) 
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Cp^  j  •  7.  01  cal/(gm  mole)(*K); 

CP,  *  C'pCH,)  +  X0,  CP«I,) 

*  (2/3)<6.88)  +  (1/3)(7.01)  *  6.9233. 
d.  Specific  Heat  Ratio 


6.9233 

6.  9233-1.98719 


1.4026. 


e.  Heat  of  Formation  6H'^|  of  Each  Reactant  at  the  Initial  Conditions 

By  definition  of  the  heat  of  formation  of  an  element 

(in  its  standard  state  of  aggregation)  at  tj  is  equal  to  zero. 

Therefore, 


A  Hi 

f(H2) 

^(02, 


AH 


«  0.0; 
•  0.0. 


Step  9.  Initial  Assumption  for  p^ 

p.,  *  20.00  atmos. 

Step  10.  Initial  Assumption  for 

T2  *  3800* K. 

Step  11.  Determination  of  Numerical  Values  of  the  Equilibrium  Constants  for  the 
Dissociation  Reactions 

From  Figures  12  through  15, 

Kn  =  •  3758  (atmos)Vl  at  T,  *  3800*K.  ,* 

KtHjO)  e- 

KL  ■  .  5984  (atmos)1^  at  T,  =  3800*K.  ,** 

PtHfo)  2 


*  This  value  also  obtainable  from  Ref.  8,  Table  4,  by  graphical  interpolation. 

**  This  value  also  given  in  Ref.  7,  page  25. 
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Kn  ■  1. 107  (atmoa)1*  at  T,  ■  3800*K.  ,* 

P(Ht>  2 


‘PtO.1 


■  1.032  (atmos)1'*  at  T,  =  3800*K.  ** 


Step  12.  Calculation  of  the  Partial  Pressures  (pj's)  in  the  Final  Mixture 
Referring  to  equations  (114)  and  (115).  let 


(116) 


;|g  -  ^KPtogl  KP<HgOl  4.  KP(Hgo)  +  A 

\  p«-  pJ  J 


(117) 


Pis  *  (PHf  +  Kp(„t,  P„  *  “  P*  >  » 


(118) 


k„  s  2 


PH- 


2  C„  , 


(119) 


x  /KP<0»>  KP<HtO)  +  KP(HgO>  \ 
\  pHt  2p  ■„*/ 


(120) 


-  (  PH,  +  3  KP(„t,  )  ' 


Then,  from  equation  (114), 


-  C(t  ±  JcJ  -  4  C|,  Cl5 


*  This  value  also  given  in  Ref.  7,  page  25. 
*•  This  value  alao  given  in  Ref.  6. 


(121) 


(122) 
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and,  from  equation  (115), 


(123) 


For  the  solution  of  equations  (12 2)  and  (123),  values  of  pH^  have  been 
assumed,  as  indicated  in  Table  1.  These  values  of  pH^  ,  in  conjunction 

with  the  values  of  the  equilibrium  constants  (Step  11)  and  the  assumed 
final  pressure  (Step  9)  have  been  used  to  determine  the  values  of  the  con¬ 
stants  cjj,  cj2»  c j 3,  andkjj,  kj^,  kjj  (Equations  (116)  through  (121)). 

Equations  (122)  and  (123)  were  then  solved  to  obtain  two  values  of  Ph20  ’ 
indicated  in  Table  1  as  p^  (Equation  (122)) and  pH2Q (Equation  (123)).  The 

difference  between  these  two  values  of  p|t^Q  is  designated  in  Table  1  as 
APh2o  *  The  values  of  ApH20  obtained  were  plotted  as  a  function  of  the 
assumed  pH^  (Figure  20).  After  two  trials  (i.  e.  ,  two  assumed  values 
of  p  ),  the  value  of  p  at  the  point  where  the  curve  crossed  the  Ap„  *0 
line,  was  taken  as  the  assumed  p^  for  the  following  trial.  The  fourth 
trial  gave  sufficiently  accurate  values  of  pH^  and  pH^0  to  satisfy  both  equa¬ 
tions  (122)  and  (123).  (An  alternate  method  of  obtaining  a  solution  of  the 
simultaneous  equations  (122)  and  (123)  is  to  plot  each  of  these  two  equa¬ 


tions  as  p^  vs  pM  ^  on  the  same  graph  and  determine  the  point  of  inter¬ 
section  of  the  two  curves.  It  may  be  found  that  in  some  cases  more  rapid 
convergence  is  obtained  by  this  method.) 
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TABLE  1 


Calculation  o£  n 

P2  -  20  atmos 
T2  .  3800°K; 


and 


from  equations  (122)  and  (123) 


K0  =  .  3758  (atmos)  ^  ;  K*0 

K(HtO)  ,  MHtO> 

*■>(,,,  =  1.107  (..mo.)  *;  Kp(#>) 


*  .  5984  (atmos)  ^  ; 
=  1.  032  (atmos)  ^  * 


Trial  1 

Trial  2 

Trial  3 

Trial  4 

PH 

Mt 

4.35 

3.00 

3.43 

3.  51 

C„ 

.0074634 

.0156918 

.012003 

.0114640 

c«« 

1.37607 

1.47476 

1.4362 

1.42989 

c.* 

-13.3412 

-15.0826 

-14. 520 

-14.4160 

^11 

.0149268 

.0313836 

.024006 

.022928 

kl2 

.232611 

.302019 

. 27463 

. 270193 

k„ 

-5.  50442 

-3.95870 

-4.4551 

-4. 54698 

PMt0(Eq-  (122)) 

9.23306 

9.  30578 

9.3768 

9.  3733 

PMfo(Eq-  (123)) 

12.9320 

7.4068 

9.0550 

9. 3732 

APHtO 

+  3.  6989 

-1.8990 

3218 

-.0001 
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Using  equation*  (106)  through  (109),  respectively,  to  calculate  the 
values  of  the  partial  pressures  of  the  remaining  four  constituents, 


^om  = 


p„ » 


/*P(H gO)  pHfO _  1 

(,3758)(9.3732)T  .  lmWnt 

\  pj,  /  ”  1 

3.51  J 

KP<h,oi  P"^  . 

P»,T 

(.5984)  (9.3732)  .  2.  „stl 

(151)'* 

KP(h,)  Ph,*  * 

1 

(1.107)  (3.51)*  *  2.07396  atrooa. 

P»,o 

(I032)(.3758)(a3732)  ,  j  0356fc  atmo. 

P", 

3.51 

Rounding  off  to  three  decimal  places  all  of  the  values  of  partial  prea- 
sure  obtained  above,  and  adjusting  the  largest  of  these  values  (  pM^0  ) 

slightly  in  order  to  give  a  total  (p^)  of  20.000,  the  following  values  of 
partial  pressure  (p.'s)  are  obtained: 

J 


\ 

S 

3.  510  atmos 

0 

M 

9.  379 

HjO 

D 

m 

1.007 

D 

s 

2.994 

oh 

* 

K 

2.074 

M 

1.036 

Step  13. 


2  pt  s  p  •  20. 000  atmos 

Total  Number  of  Moles  of  the  Final  Mixture 
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Using  equation  (75), 


n2  =  P2 


Jn'iMw(i) 
1 _ 

2pj  MW(j) 


=  p. 


nHg 


+  n0*  M#(oa> 


PH^>  M«<«,0l  +  P„,  “»(»,!  +  +  P0„  “"tOHl  +  Ph  M»(Hl  +  P0  M»(0) 


20.000 


(I)  (2.0160)  +  (.5)  (32.0000) 


(».3T6)(ieOI«0)  +  (3.9IO)(2.0l60)+(l.007)(32.0000)-H^994)(l7.0080H-(2.0?4)(L0080)+(l.036)(l6.000) 


or 


n  «  1. 29676 

w 


Step  14.  Composition  of  the  Final  Mixture 


Using  equation  (73)  and  values  of  p.  calculated  in  Step  12, 


(9.379)  (1.29676) 

20.000 


.  608116 


(3.5)0)0.29676) 

20.000 


. 227581 


n0* 


(1.007)0.29676) 

- 551*00 — 


.065292 


* 


(2.994)  (1.29676) 
20000 


n0 


(2.074)  (1.29676) 

20.000 


(1.036)0.29676) 

•  - 26.5W" 

fni  *  "2 


«  .194125 

s  .134474 

•  .067172 

*  1.29676 
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Step  1 5.  Propertiea  of  the  Final  Mixture 

a.  Mole  Fractions 


*HjO 

m 

sr 

3 

kh2o 

X 

_  9.379 
'  20.000 

• 

. 46895 

Xu 

nH2 

■ 

*82 

3  3,310 

* 

. 17550 

n* 

"2 

20000 

*0. 

5 

n2 

= 

fT 

1.007 

20.000 

« 

.05035 

Xoh 

n0H 

m 

m 

.  2.994 

= 

. 14970 

n2 

P2 

20.000 

Xh 

nH 

>JL 

r  ?,07±_ 

r 

. 10370 

n2 

^2 

20.000 

X0 

s 

no 

n2 

m 

p. 

.  1.036 

20.000 

= 

.05180 

= 

1.00000 

j 

b.  Molecular  Weight 

M"Z  *  ^  Xj  M*U> 

=  Xh2o  MW(|fg0)  +  Xh2  Mw(M2)  +  Xo2  MW(02, 

+  Xoh  +  X„  MW(h^  +  X0  Mwtoi 

=  (.  46895)(18.  0160)  +  (.  17550)(2.  0160) 

+  (.  05035)(32.  0000)  +  (.  14970)(17.  0080) 
'  +  (.  10370)(1. 0080) +  (.  051 80) (16.  0000), 

or 

MW2  »  13.  8930. 

c.  Molal  Heat  Capacity* 

At  T2  =  3800®K, 

Cp^oj  =  13.  87  cal/(gm  mole)(°K), 

Cp(H  ^  «  9.  10  cal/ (gm  mole)(°K), 

Cp^o2)  *  9.  85  cal/(gm  mole)(°K), 

*  The  molal  heat  capacities  for  H2O,  H2,  O2  and  OH  were  obtained  from  Figure  16 
through  19,  respectively;  and  the  values  for  O  and  H  were  calculated  from  enthalpy 
data  given  in  Ref.  7,  page  22. 
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cp<»o 

a 

9.115  cal/(gm  mole)(°Kj, 

CP  tHt 

s 

4.965  cal/(gm  mole)(*K), 

a 

4.965  cal/(gm  mole)(°K). 

CPt  =  ^  "i 

CP|J» 

a  XM*0 

CP<»^  +  X«. 

CP(H2»  +  X02  CP(0jj)  +  X0H  CP(OH) 

+  XH 

CP.H,  +  *■> 

CPW 

•  {.  -46895) (13.  87)  +  {.  17550)(9.  10)  +  (.  05035)(9.  85) 

+  (.  14970)(9.  115)  +  (.  10370)(4.965)  +  (.  05180)(4.  965). 


Cp2  s  10.  7339  cal/(gm  mole)(*K). 
d.  Specific  Heat  Ratio 


10.7339 

10.  7339-1.98719 


1.  2272. 


e. 


Heat  of  Formation  AH 


of  Each  Product  at  the  Initial  Temperature1*1 


By  definition  of  AH'^,  in  the  same  manner  as  in  Step  8  for  the 

initial  mixture,  the  heat  of  formation  of  an  element  (in  its  standard 
state  of  aggregation)  at  is  equal  to  zero.  Therefore, 


AH',  »  o.o, 


A  H* 


f<ot»  * 


0.0. 


From  the  literature,  at  t_ 


AH,  ■  +57,801** 

' L  (HgO) 


:  18°C, 
cal/gm  mole, 


but, 


since  t 

L 


tj  in  this  case, 


A  Hi  =  AHf,  . 

MH^S)  *  t*  (HjO) 

*  See  footnote  in  Step  8,  par.  e. 

**  Ref.  2,  page  255. 


*  +57,801  cal/gm  mole. 
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From  the  literature,  at  t^  •  25*C, 


AH 


fL  (OM)  *  -9,  310*  cal/gm  mole, 


*»  -51, 900*  cal/gm  mole, 


AH 


*  -59, 100*  cal/gm  mole, 


but,  since  t  wj*  t  in  this  case,  Equation  (81)  is  used  to  calculate 
L  1 

values  of  AH*  for  OH,  H  and  O**: 

Mj) 


AH'<o>  ■  A%.  -H'i+  • 


(81) 


Values  of  H'j  and  £H'e...  obtained  from  the  literature  for  OH,  H 


SU) 

and  O  are  indicated  in  Table  2. 


Frost  the  shore  data  and  equation  (81), 


AH'*«OH»  -  +2H‘. 


(OH) 


r  -9,  310  +  47.  9  -  48.  5  «  -9310.  6  cal/gm  mole; 


ah; 


f(H) 


aVh,-h;h.  + 


(H) 


a  -51 , 900  +  33.  2  -  24.  0  =  -51, 890.  8  cal/gm  mole; 


AH 


f(0> 


AHtu.»  -H«„  +2h’.<.. 


b  -59,  100  +  33.  2  -  24.  5  «  -59,  091.  3  cal/gm  mole. 


*  Ref.  7,  page  26. 

**  The  difference  between  AH’f^  and  AHf^j  for  this  case  (t^  -  tj  =  7°C)  is  negli¬ 
gible  (<0. 1%).  The  calculation  of  A  H'f^  according  to  equation  (81)  has  been  included 
for  purpose  of  illustration. 
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TABLE  2 


Literature  Values  of  Enthalpies,  cal/gm  mole  of  "j" 


Component 

"j" 

Hi 

(cal/gm  mole) 

^H'«(j) 

(cal/gm  mole  of  "j") 

Reference 

OH 

-47.  9 

-48.  5 

Calculated  (*)  from 

H 

-33.2 

-24.0 

data  in  Ref.  7, 

0 

-33.  2 

-24.  5 

page  22  and  Ref.  2, 

page  216. 

*  The  calculation  of  H' •  and  2hL...  from  data  obtained  in  References  2  and  7  is  shown 

J  e0) 

below  for  OH. 


H- 8  /©,«, dT=  A-  dT  -  (V\), 


(124) 


where  Cp 
therefore,  H' 


cal 


■  mean  Cn  between  Tt  (298° K)  and  T,  (291° K)  *  6.  8421  7 - , 

(OH)  P  L  1  (gm  mole)(°K); 


tOH)  *  (6.  8421)(291-298)  a  -47.9  cal/gm  mole. 

The  formation  reaction  for  one  mole  of  OH  is  1/2  H £  +1/2  O2 . -»  OH;  then, 

^*H«(OH)snH2  £P(h2)  (V\)  +  no8  CPf0a)  (T.-TJ  , 

where  Cp^j  •  mean  Cp  between  T^  (298°K)  and  Tj  (291*K)  *  6.86  cal/(gm  mole)(*K), 
and  CP{0  j  *  mean  Cp  between  T^  (298° K)  and  Tj  (291*K)  =  7.00  cal/(gm  mole)(0,K); 

2  H*  =  (l/2)(6.86)(291-298)+l/2(7.00)(291-298)  =  -48.  5  ~ - - - — -  • 

"  8  (oh)  gm  mole  of  OH 

The  values  of  the  heats  of  formation,  AH'f^j,  calculated  above,  are  summarized 
for  convenience  in  Table  3. 
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TABLE  3 


Heats  of  Formation  of  Products,  cal/gm  mole 


AHf(H20) 

0.0 

0.0 

+  57,801 

-9310.6 

-51,890.8 

-59,091.  3 

Step  16.  Heat  of  Reaction,  Energy  Release  per  Mole  of  Initial  Mixture 
Using  data  from  Step  14  and  Step  15  (par.  e), 

f  ni  n»KAHw«+n".AHw+  n°«AHw+  n»»AHw+  nM,H,  +n°  AH»  to. 

*  (.  6081l6)(+57,  801)  +  (.  227581)(0.0)  +  (.065292)(0.0) 

+  (.  1941 25)(-9310.  6)  +  (.  1 3447 4) (-  51,  890.  8)  +  (.  067172) 
(-59,091.  3) 

=  +22,  395.0  calories; 


n, 


*  nM2  AHf{M2)  +  n®2  AHf(02) 

*  (1)  (0.0)  +  (1/2)  (0.0)  =  0; 

=  n'Hz  +  :  1+2  *  1.5  moles  % 


then,  by  equation  (83), 


?”j  AHi(i)  •  ?">  AH’f<i 


(i) 


AH, 


n. 


+  22,  395,0  -  0.0 
1.5 


or  AHj^  s  +14,  930.0  cal/gm  mole  of  initial  mixture. 

Step  17.  Energy  Release  Function 

Using  equation  (84)  and  data  from  Steps  7,  8  (par.  d)  and  16, 
AHR  _  +14, 930.0 _ 


B  = 


Y.RT, 


( 1 . 40 26) ( 1 . 9871  9)(291) 


=  18.407; 
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Step  18,  Gamma  Ratio 


Step  19. 


Step  20. 


Step  21. 


Step  22. 


Step  23. 


Step  24. 


From  Figure  5  and  data  from  Steps  8  (par.  d.)  and  15  (par.  d. ) 
or  by  direct  calculation. 


1.2272 
*  1.4026 


.87495  . 


Molecular  Weight  Ratio 


Using  data  from  Steps  8  (par.  b  )  and  15  (par.  b  ), 


Mw2 


13.  8930 
12.0107 


1.15672  . 


The  Dimensionless  Parameter 


E 


From  Figure  1  and  data  from  Steps  15  (par.  d  )  and  17,  or  using 
equation  (47), 

E*  (Y*-|)B  =[(1.  2272)2-l]  [jl8.407]  =  9.  314'31  . 

Gamma  Function 


From  Figure  2  and  data  from  Steps  8  (par.  d  )  and  15  (par.  d  ), 
or  using  equation  (48) , 

r  =  — ~ — 3 —  s  (!»  ZZJZ)Z_-  1^4026  _  0.  183146  . 

(KJ  -|)  1.4026(1.4026-1) 


Pseudo-Mach  Number  (Squared) 

From  Figure  3  and.  data  from  Steps  20  and  21,  or  using  equation 

(45), 

M*  s  2  E  +  2f  ■  (2)(9.  31431)42(0. 183148)  =  18.  99491  . 


of  M|^  and  ^ 


from  Steps  18  and  22, 


Mach  Number  (Squared) 

From  Figure  4,  and  values 

Mj^  *  19.0  (Same  as  Rd^  within  the  accuracy  of  reading  Figure  4). 

The  Dimensionless  Parameter  T,  M 

From  Figure  9  and  data  from  Steps  8  (par.  d  )  and  23,  or  by 
direct  calculation, 
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T,  Mj2  •  (1.4026)09.0)  ■  26.65  . 

Step  25.  Temperature  -  Molecular  Weight  Ratio 

From  Figure  8  and  data  from  Step*  15  (par.  d  )  and  24,  or  using 
equation  (85), 


T,  M. 


Step  26.  Final  Temperature 

Using  data  from  Steps  7,  19  and  25, 

T2mwi  mw? 

T1  am  ■  •  -  -  «  T 

2  TjMw.,  A1 

=  (7.0972)(1.  15672)(291) 

*  2389° K . 

Since  this  calculated  value  of  did  not  check  the  assumed  value 

(Step  10),  other  values  of  were  assumed  and  Steps  10  through  26 

were  repeated  until  sufficient  data  were  obtained  to  permit  graphical 
determination  of  the  correct  value  of  T^  corresponding  to  the  assumed 
pressure. 

k  plot  of  A  PjjgO  TS  p^  (uauatd)  for  sash  assuasd  T2  la  shorn 

In  figure  21,  and  a  plot  of  T2  (9*1^^)  r»  T2  (>MUM>g)  1*  chova  la 

figure  22  (for  Trials  1  through  29,  l.e»,  for  pg  ''rimml)  *  20»00 

•taos).  Table  4  Hats  prlnoipel  lntenasdiate  ralues  oaleulated  la  those 
trials. 

The  eorrect  value  of  T2  corresponding  to  the  first  assaned  ot 

20.00  ataos  Is  shovn  In  figure  22  to  be  35^5 (tbs  point  of  later* 
section  of  the  curve  vlth  the  45°  lias).  This  value  of  t2  Is  also  Indi¬ 
cated  la  Table  4  (Trials  25-29)  as  T2  (adjusted)* 


(Y*+i)«  ("&  M*  +  X  «i*  +  2) 

1.  2272  .  1 

1  1  (26.  65  +  /  c  +  2) 

(1.  2272+I)2  26,65 


7.0972  . 
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It  was  considered  that  the  value  of 


^(adjusted) 


obtained  in  Trial 


29  as  described  above  was  so  close  to  the  value  of 


T 


2(as  Burned) 


that 


further  calculations  were  not  warranted.  This  was  considered  to  be 
justifiable  since  the  thermodynamic  data  (K^  and  C^)available  were 

not  of  sufficient  accuracy  to  give  data  different  from  that  obtained  in 
the  last  trial  (Trial  29).  Hence,  all  intermediate  values  calculated  in 
the  last  trial  (Trial  29),  based  on  the  assumed  T 2  for  that  trial,  for 
all  practical  purposes  can  be  considered  to  be  accurate  values  corres 
ponding  to  the  adjusted  T2. 


Step  27. 


Pressure  Ratio 


Using  the  values  of  and 
ure  6  or  using  equation  (87), 


from  Trial  29,  then  from  Fig- 


s  Ta'+i  + 

*  l'.'aml+i  <”• 3793  +  » 

=  17.  3887  . 


Step  28.  Final  Pressure 

Using  data  from  Steps  6  and  27, 


p2  =  —  •  pt  ■  (17.  3887)(1)  «  17.3887  atmos. 


Since  this  calculated  value  of  p^  did  not  check  the  assumed  value 

(Step  9),  other  values  of  p  were  assumed  and  Steps  9  through  28  were 

repeated  until  sufficient  data  were  obtained  to  permit  graphical  deter¬ 
mination  of  the  correct  value  of  p  . 


A  plot  of  ApHz0  vs  PH2<a8BUmed)  for  each  a8Sumed  T2  is  shown 

in  Figures  21,  23,  25,  26  for  assumed  p^'s  of  20,  17,  17.  17  and  17.  20, 
respectively.  A  plot  of  T2(calculated)  vs  T2(a88ume<J)  is  shown  in 

Figures  22,  24  and  27  for  assumed  P2's  of  20,  17  and  17.  20,  respectively*. 


*  For  an  assumed  p2  of  17. 17,  no  plot  of  T2<calculated)  vs  T2<a8Buined)  was  necessary 
since  the  calculated  T2  for  the  first  trial  was  so  close  to  the  assumed  T2  that  additional 
trials  were  not  warranted. 
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Step  29. 


Step  30. 


Step  31. 


Figure  28  .hew.  .  plot  of  P2(calculate<1)  v.  P2(„<u  ned).  Table  4 
lists  principal  intermediate  valuea  calculated  in  Trials  1  through  53. 

The  correct  value  of  the  final  detonation  pressure,  p  ,  is  shown 

2 

in  Figure  28  to  be  17.  22  atmos  (the  point  of  intersection  of  the  curve 

with  the  45*  line).  This  value  of  p.  is  also  indicated  as  p„.  . 

2  r  2  (adjusted) 

in  Table  4. 

It  was  considered  that  the  value  of  P2(adjugted)  °kta*ne<*  *n  Trial 

53  as  described  above  was  so  close  to  the  value  of  p_.  that 

r2(assumed) 

further  calculations  were  not  warranted.  Hence,  all  intermediate 
values  calculated  in  the  last  trial  (Trial  53),  based  on  the  assumed  p^ 
for  that  trial,  can  be  considered  to  be  accurate  values  corresponding 
to  the  adjusted  p^. 

Detonative  Mach  Number 

Using  the  value  of  M^  from  Trial  53, 

Ml  *  -  ^2 6.  40  =  5.  13809  . 

Detonation  Velocity 

Using  equation  (89)  and  data  from  Steps  7,  8  (pars,  b  and  d  ) 
and  29, 


D  *  M.  A 7'  «  T, 

1  V  MW| 


=  5.  13809 


/( 4.  185x  107)(1.4026)(1.98719)(291) 
V  12.0107 


■  2.73149xl03  cm/sec. 

This  detonation  velocity  expressed  in  more  common  units  is 

2731.  5  meters/sec 


or 


8961.  6  ft/sec 

Calculation  of  Additional  Data 
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a.  Density  Ratio 


E 


P  ?  T  iM*. 

Using  the  value*  of  —  and-- - ±-  from  Trial  53, 

PX  TjH- 


then  from  Figure  10  or  using  equation  (15), 


Pz  P2  /  T2Mw1  17.2269 
pt  *  /  T1Mw2  “  9.67764  *  K  78007 

b.  Initial  Density 

Using  equation  (90)  and  the  values  of  P, ,  T  and  from 

11 

Steps  6,  7  and  8  (par.  b  )  of  this  sample  calculation, 

„  MwlPl  (12. 0107) (1.0133x10^)  _  _ 4  ,  , 

°  *  JRT j  *  (4. 185 xl07)(l.  98719)(291)  *  5-02895xl0  gm/cm3 

or  initial  density  is  0.5029  gm/liter. 

Final  Density 

A 


c. 


Using  the  above  values  of  *j3.  and  p 


P* 


•  f>  *  (1.78007)(5.02895xl0-4) 

n 

•  8.95188x  10**  gm/cm3, 
or  final  density  is  0.  8952  gm/liter. 
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Density  Ratio 


P>  TjMwi 

Using  the  values  of  —  and-——- — ±-  from  Trial  53, 

Pj  TjMw^ 

then  from  Figure  10  or  using  equation  (15), 


17. 2269 
9.67764 


1.78007  . 


Initial  Density 


Using  equation  (90)  and  the  values  of  P^,  and  from 
Steps  6,  7  and  8  (par.  b  )  of  this  sample  calculation, 


ft 


MwlPl  (1 2. 0107)(1. 01 33  x  10^)  . ^4  , 

*  JRTj  *  (4.185x107){1.98719)(291)  *  5- 02895x10 


or  initial  density  is  0.5029  gm/liter. 


Final  Density 


A 


Using  the  above  values  of  and  p  , 

“\  1 

f>  z  A.  .  f>  •  (1.78007)(5. 02895  xlO'4) 

P\ 

m  8.  95188xl0"4  gm/cm3, 


or  final  density  is  0.8952  gm/liter. 
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FIGURE  25.  DETERMINATION  OF  COMPOSITION  OF  FINAL  MIXTURE,  TRIALS  45  THROUGH  49. 
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FIGURE  26.  DETERMINATION  OF  COMPOSITION  OF  FINAL  MIXTURE,  TRIALS  50  THROUGH  53. 
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